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Higher orders of the high-temperature expansion for the Ising model in 
three dimensions 
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The new algorithm of the finite lattice method is applied to generate the high-temperature expansion series of 
the simple cubic Ising model to /3 for the free energy, to /3 for the magnetic susceptibility and to P for the 
second moment correlation length. The series are analyzed to give the precise value of the critical point and the 
critical exponents of the model. 



1. INTRODUCTION 



The finite lattice method '1'2"3' is a power- 
ful tool to generate the exact high- and low- 
temperature series and other exact expansion se- 
ries for the spin models in the infinite volume 
limit. Using this method we can skip the job of 
listing up all the relevant diagrams and of count- 
ing the number they appear, which is inevitable in 
the graphical method. Then the main task is re- 
duced to the calculation of the partition function 
for the appropriate finite size lattices, which can 
be done efficiently by the site-by-site integration 
using the transfer matrix formalism |4I5| without 
the graphical technique. The CPU time and the 
computer memory needed to obtain the series to 
order N increase exponentially with N in two di- 
mensions, while they grow up exponentially with 
A^^ in three dimensions. Thus the finite lattice 
method has been extremely effective in two di- 
mensions, but it has generated the series only in 
some limited cases in three dimensions. 

Recently we presented the new algorithm of 
the finite lattice method '6l7j in which the CPU 
time and the computer memory increase expo- 
nentially with A^logiV. It enabled us to extend 
the high-temperature expansion series for the free 



energy of the simple cubic Ising model to order 
/3^^ from the previous series of order which 
was given by the previous algorithm of the finite 
lattice method, and to order (3^"^ for the magnetic 
susceptibility from the previous series of order 
/?^^, which was given by the graphical method. 

Here we apply the method newly to the high- 
temperature expansion of the second moment cor- 
relation length to order /?^^ and also extend the 
series for the free energy to order /3^°. The series 
themselves will be presented elsewhere and we 
give here the result of the analysis of the newly 
generated series combined with the series for the 
magnetic susceptibility already obtained. 

2. SUSCEPTIBILITY AND CORRELA- 
TION LENGTH 

The magnetic susceptibility x and the second 
moment correlation length ^2 are expected to be- 
have like 
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6' - A2 {(3c ~ (3)-'" [1 + C2{Pc - f3)^ + ■ ■ ■ ] 

Here /3c is the critical point and A is the leading 
correction-to-scaling exponent. 

There are several evaluation of the value of A to 
be around 0.5. To make more precise evaluation 
we plot 2z/ versus A in Figure Q in the analysis of 
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Figure 1. Critical exponent 2v plotted versus A 
in the M2 method of the series analysis for £_2^ . 



the ^2 series by the M2 inethod|HI- The narrow 
crossing region of the lines gives A = 0.51(1). 
The modified ratio series [TD] 
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made of the expantion coefficients a„ for each 
quantity is expected to behave like 



(/3c)„ = /3c + 



,1+A 



(2) 



Although the modefied ratio series for the longest 
series of each of the x and can be fitted well 
by the first two terms in Eq. (O with A = 0.51, 
there is still a considerably large ambiguity in ex- 
trapolating the plot to n —>^ 00. 

Then we take the following combined quantity 
of the susceptibility and the correlation length. 



.(3) 



We can expect that, if we take an appropriate 
value of 07, the third and later terms are very 



Figure 2. Modified ratio series for the combina- 
tion of X E^nd ^2 with uj — 0.055 plotted versus 
l/ni+^ 



small and the ratio series Q for this quantity can 
be fitted quite well by the first two terms in Eq. 

In fact we find as in Figure |21 that the situ- 
ation is realized for co = 0.055 for even series of 
n = 22-28 giving Pc = 0.2216545(1) if we assume 
A — 0.51(1). Similarly the combined quantity 
6Vx"' with uj' = -0.15 gives /3c = 0.2216542(4) 
with the same value of A assumed. 

There are two estimates for the critical point 
of the model by the recent large scale numeri- 
cal simulations. One is /3c = 0.22165459(10) by 
Blote et a/. 22 using the cluster algorithm and 
the other is /3c = 0.2216545(15)11' by Ito et 
al. using the non-equilibrium relaxation method. 
They are contradictory to each other. Our esti- 
mation by the analysis of the high-temperature 
series is consistent with these estimates. 

As for the critical exponents the modified ra- 
tio series for the combination (PJ of the magnetic 
susceptibility and the second moment correlation 
length gives 7 = 1.2369(3) and v = 0.6298(3) for 
A = 0.51(1). These values are consistent with 
the estimate of 7 = 1.2373(2) and v = 0.6301(2) 
by the high-temperature expansion of the gener- 
alized Ising model J3 ^'^^ of 7 = 1.2396(13) and 
V — 0.6304(13) by the perturbative renormaliza- 
tion method[n|. 
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3. SPECIFIC HEAT 

For the critical exponent of the specific heat 



= A4i(3c-Pr''[l + c4{f3c-f3)^ + --- ], 
the ratio series 

0-271-2 



with n < 50 gives a = 0.1035 ± 0.0005 for 
f3c = 0.2216546(1) and A = 0.51(1). This esti- 
mate of a is consistent with that by the inhomo- 
geneous differential approximation for the high- 
temperature series. From the hyperscaling rela- 
tion a = 2 — du, the range of = 0.6295 — 0.6305, 
within which all of the recent precise estimates 
of the critical exponent for the correlation length 
falls, gives a = 0.1085 - 0.1115. Thus there is a 
discrepancy between the estimate of a from the 
high-temperature series of the specific heat and 
the estimate from the hyperscaling relation. 

4. SUMMARY AND DISCUSSION 

We have calculated the high-temperature se- 
ries for the second moment correlation length of 
the simple cubic Ising model using the new al- 
gorithm of the finite lattice method. Combining 
this series with the previously obtained series for 
the magnetic susceptibility, we have given pre- 
cise estimates of the critical point and the critical 
exponents, which is consistent with other recent 
estimates. We have also extended the series for 
the free energy, which gives the estimate of the 
critical exponent for the specific heat, which is 
inconsistent with that derived from the critical 
exponent of the correlation length using the hy- 
perscaling relation. 

The new algorithm can also be applied to the 
high-temperature expansion of the models that 
have continuous spin variables such as the XY 
model in three dimensions. Especially it is im- 
portant to calculate the high order series of the 
high-temperature expansion for the XY model, 



since the phase transition of this model belongs 
to the same universality class as the lambda tran- 
sition of the helium and the specific heat exponent 
at this phase transition point was observed very 
precisely in the Space Shuttle experiment 15 . 

Furthermore the basic idea of the new algo- 
rithm can be used in the low-temperature expan- 
sion for the spin models. We can expect that the 
new algorithm will enable us to generate the se- 
ries for these models that are much longer than 
the presently available series. 
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